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ABSTRACT 

In order to understand the physical mechanisms underlying non-steady stellar spiral arms in disk 
galaxies, we analyzed the growing and damping phases of their spiral arms using three-dimensional 
A-body simulations. We confirmed that the spiral arms are formed due to a swing amplification 
mechanism that reinforces density enhancement as a seeded wake. In the damping phase, the Coriolis 
force exerted on a portion of the arm surpasses the gravitational force that acts to shrink the portion. 
Consequently, the stars in the portion escape from the arm, and subsequently they form a new arm 
at a different location. The time-dependent nature of the spiral arms are originated in the continual 
repetition of this non-linear phenomenon. Since a spiral arm does not rigidly rotate, but follows the 
galactic differential rotation, the stars in the arm rotate at almost the same rate as the arm. In other 
words, every single position in the arm can be regarded as the co-rotation point. Due to interaction 
with their host arms, the energy and angular momentum of the stars change, thereby causing the 
radial migration of the stars. During this process, the kinetic energy of random motion (random 
energy) of the stars does not significantly increase, and the disk remains dynamically cold. Owing 
to this low degree of disk heating, the short-lived spiral arms can recurrently develop over many 
rotational periods. The resultant structure of the spiral arms in the TV-body simulations is consistent 
with some observational nature of spiral galaxies. We conclude that the formation and structure of 
spiral arms in isolated disk galaxies can be reasonably understood by non-linear interactions between 
a spiral arm and its constituent stars. 

Subject headings: galaxies: kinematics and dynamics - galaxies: spiral - galaxies: structure - methods: 
numerical 



1. INTRODUCTION 

The spiral arms of disk galaxies are the most promi- 
nent structures, and the arms are formed due to gravita- 
tionally driv en variations in the surface density i n the 
stellar disk flRix fc Zaritskvl 119951: iGrosbol et ap l2QQ4|; 
IZibetti et all 120091: lElmegreen et all 1201 If ) . It is known 
that the spiral arms of disk galaxies can be excited by 
tidal interactions with nearby companion galaxies ( e . g . , 
lOh et al.l[2008l:lDobbs et aTJl20int IStruck et al.ll201lD. as 
well a s by the central stellar bar (e.g.. lSellwood fc Sparkd 
119881 ). In addition, spiral arms can also be self- induced 
and maintained without these external gravitational per- 
turbations, and further, they can propagate stationary 
de nsity waves in glo bally stable disks, as hypothesized 
by iLin fc Shul (|1964[) . However, the physical origin and 
dynamical evolution of spiral arms in disk galaxies have 
thus far not been fully understood. In this study, we 
focus on the dynamics of stellar spiral arms and stars 
around them in a disk galaxy without external perturba- 
tions and a bar structure. 

The Lin-Shu theory posits that the stellar spiral 
arms can be interpreted as (quasi-) s tationary density 
waves with a cons tant pattern speed (|Lin fc Shul 119641 : 
iBertin fc Lin|[l996l ): the theory has gained wide accep- 
tance as providing a clear explanation regarding the dy- 
namics o f spiral arms. However, as pointed out in the 
study by iToomrel (|l969f ). the quasi-stationary hypothe- 
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sis has a serious limitation because of a dispersal nature 
of the tight winding spiral waves, they radially propagate 
with the group velocity, and are consequently absorbed 
at t he Lindblad resonances (|Lvnden-Bell fc Kalnaisl 
1972). Therefore, the stationary waves that last for 
a long duration re quire some a mplification mechanisms 
such as WASER (|Markj 119761) and a feedback cycle 
that involves the reflection of the inward propagating 
wave into an outwa rd propagating one at the Q-barrier 
(jBertin et al.lll989aH bT). 

However, nearly all the previous time-dependent sim- 
ulations that have been executed thus far have been un- 
able to prove the existence of stationary density waves in 
a disk galaxy without external per t urbat ions and a bar 
structure 3 . ISellwood fc Carlberel (|1984D claimed that 
self-induced stellar spiral structures in stellar disks arc 
not stationary and fade away within an interval of about 
ten rotational periods. Their study stresses on the signif- 
icance of cooling mechanisms, such as dissipation by the 
interstellar medium (ISM), that are required to main- 
tain stellar spiral arms. More recently, iFuiii et al.l (|2011f ) 
used 3-D A-body simulations to show that although stel- 
lar spiral arms are short-lived, they are also recurrently 
formed, and as a result , the spiral featu res are maintained 
over 10 Gyr (see also iSellwoodl l201lT ) . This is because 
of negative feedback that causes the dynamical heating 
of stars in the disk galaxies due to spirals. The non- 

3 IThomasson et al.l lfl990l 'l and lElmegreen k. Thomassonl lTl993h 
have suggested that the feedback cycle caused by the reflection at 
the Q-barrier causes the long life of the m = 2 — 3 spiral arms. 
However, their ./V-body spirals are not stationary; instead, a timc- 
dependen t evol ution of the spiral features is observed. See also 
ISellwoodl |(20TJ). 
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steady spirals are also subject to a similar phenomenon 
when the dynamics of the ISM are self-consistently solved 
with stars in disk galaxies. Both the stellar spirals and 
the ISM undergo mo tion in the wake of the the galac- 
tic rotation (see also iGrand et al.l l2012d ). In fact, the 
spirals can be considered as being "wound". The ISM 
forms dense regions associated with the non-steady stel- 
lar spirals; however, th ese regions are no t the conven- 
tional "galactic shocks" ([Wadaet al.l 120111) . 

The existence of non-steady stellar spirals has been 
suggested by certain observations in our Galaxy. For 
example, the age-velocity di spersion relation of st ars in 
the solar neighborhood (e.g.. lHolmberg et al.ll2007f ) most 
naturally accounts for the exi s tence o f the non-steady 
spira l s (ICarlberg fc Sellwood) 119851: iBinnev fc Lacevl 
19881: Uenkins fc Binnevl Il990t iDe Simone et alj I2004D . 
Baba et al.l ( 20091 ) analyzed the kinematics of star- 



forming regions using numerical simulatio n data of a 
barre d galaxy such as the Milky Way galaxy (jBaba et alj 
120101 ) . and they found that the kinematics of these re- 
gions are consistent with the observed peculiar motions 
of maser sources. These studies support the view that 
the Galactic stellar spiral arms are non-steady instead of 
stationary density waves. 

However, the current theoretical understanding of the 
dynamics of non-steady spirals is thus far insufficient. It 
has been suggested that a key factor towards understand- 
ing their dynamics is the orbital evolution of stars. In 
this study, we analyze the orbits of stars associated with 
non-steady spirals in the growth and damping phases. 

In Section [3J we summarize the numerical model 
and the me t hod, which are the same as reported by 
IWada et al.l (j2011[ ). except that the present model does 
not include the ISM, i.e., our study involves pure N- 
body simulations. We examine the global evolution of 
the stellar disk in Section (3[ in which we show the non- 
steady nature of spirals developed in the galactic disk. In 
Section 2J we focus on a typical spiral arm and describe 
its dynamical evolution, i.e., the amplification (Section 
14. ip and destruction (Section 14. 2[) processes. Section [5] 
examines the evolution of stars around the growing spi- 
ral arm from the viewpoint of the angular momentum- 
energy space. Finally, in Section [6j we compare our sim- 
ulation results with observations, and we provide a per- 
spective view on grand-design spirals. The grand-design 
spirals (i.e., m = 2 spirals) 4 will be the subject of our 
future studies. 

2. MODELS AND METHODS 

The simulations of the formation and evolution of spi- 
ral arms were carried out using a three-dimensional TV- 
body simulation. The initial point of our simulation was 
a disk galaxy in a nearly equilibrium state. The density 
profile of the stellar disk in this case is given by 

p sd (R, z) = ^ Isd 2 cxp(-R/R sd )sech 2 (z/z sd ), (1) 

47T# sd z sd 

where M s d is the mass of the stellar disk, i? sd is the radial 
scale-length, and z sd is the vertical scale-length. Since 

4 Although the term grand-design spiral indicates a reasonably 
cohe rent and extensive spiral arm in the stellar mass distribution 
(Sec Elmcgrccn & Elmcgrccn 1982, for a more precise definition), 
in a majority of cases, this means that the galaxy has m = 2 spiral 
arms JKendall et alj|20TJ) . 



the distribution function of a stellar disk is unknown, 
we generated the equilibrium state of the stellar disk by 
using an empirical method: We first generated a stellar 
disk in a near-e quilibrium state based on a Maxwellian 
approximation (jHernquistl I1993T ) , and we subsequently 
allowed the disk t o evolve for 6 Gyr under th e constraint 
of ax isymmetry (|McMillan fc Dehnenl l2007t iFuiii et all 

unni). 

We modeled a fixed potential of the dark matter 
DM) halo whose den sity profile follows the NFW profile 
Navarro et!d1ll997t) : 



Po 



where 



Po 



r/r s (l + r/r s ) 2 ' 
Cnfw 3 



4nR h 3 ln(l + C NF w) + C N fw/(1 + Cnfw) : 



C, 



NFW 



= Rh/r 



(2) 

(3) 
(4) 



an the terms Cnfw> -^h, and i?h denote the halo concen- 
tration parameter of the halo profile, the DM halo mass, 
and the virial radius, respectively. The model parameters 
adopted here are summarized in Table [T] The circular 
velocity, velocity dispersions, Q-values and r-profile of 
the stellar disk are shown in Figure [1] Here, V denotes 
a dimensionless shear rate of the disk as defined by the 
expression — dlntl/dlnR, where £1 denotes the angular 

frequency. 

We used the simulation code ASURA ([Saitoh et al.ll2008l 
120091) . where the self-gravity of stellar particle s is cal- 
culat ed by the Tree with the GRAPE method (|Makinol 
1199 1|) . We used a software emulator o f GRAPE, the 
Phantom-GRAPE dTanikawa et al.ll2013D . The opening 
angle was set to 0.5 with the center-of-mass approxima- 
tion. To perform a time integration, we used a leapfrog 
integrator with variable and individual time-steps. We 
used a Plummer softening length of e = 30 pc. This 
value is sufficiently small to re solve the three-dimensional 
struc ture of a disk galaxy (|Hernquistl 119871 : iBottemal 
l2003j) . It is to be noted that the formation of the spiral 
arms is inhibited if the adopted softening length is con- 
siderably larger than the disk thickness (~ 300 pc). The 
number of particles is 300 million. 



TABLE 1 

Model parameters for each mass component in TV-body 
simulation 



Component 


Parameter 


Value 


Dark halo 


Mass (M h ) 


6.3 


x 10 11 M Q 


(Rigid) 


Radius (i?h) 


122 


kpc 




Concentration (Cnfw) 


5.0 




Initial stellar disk 


Mass (M Bd ) 


3.2 


x 10 10 M 


(Live) 


Scale length (-R s< j) 


1.3 


kpc 




Scale height (z s( j) 


0.3 


kpc 



3. GLOBAL EVOLUTION OF SPIRALS 

Figure [5] shows the time evolution of the stellar disk. 
Hereafter, we use a rotational period measured at R = 
2i? sd (= 8.6 kpc), T rot , which corresponds to 325 Myr. 
The top panels show the face-on views of the stellar disk. 



Dynamics of Non-Steady Spiral Arms in Disk Galaxies 3 




5 10 15 5 10 15 5 10 15 

R [kpc] R [kpc] R [kpc] 



Fig. 1. — (a) Initial circular velocity curves (solid line: circular rotation, dashed: line dark matter, dot-dashed line: stars), (b) Stellar 
velocity dispersions (solid line: ctr, dashed line: ct^,, dot-dashed line: cr z ). (c) Initial Q- and T-profiles. The term T denotes a dimensionless 
shear rate of the disk as defined by the expression — dlnQ/dlnR, where Q denotes the angular frequency. The dashed curve presents the 
final Q-profile (T rot = 15). 



The middle panels show the radial distributions of the 
spiral modes which are analyzed by performing a one- 
dimensional Fourier decomposition of the disk surface 
density using the polar coordinates (R,(j)): 



n(R, 



MR) 



= J2 A m{R) C0S{77# - 4> m (R)}}. (5) 



Here, A m and (j) m {R) denote the Fourier amplitude 
and phase angle for the m-th mode, respectively (e.g, 
iRix fc Zari tskv 1991J). Upon execution of the simulation, 
spiral patterns initially developed from the noise with an 
e- folding time of about 4 galactic rotations (T rot ~ 4), 
and subsequently, the patterns settled to a nearly con- 
stant level, consistent with the report by iFujii et all 
(|201ll ). While the global features were in the quasi- 
steady state, the local spiral features were not static. 
The most prominent mode always changed on a rota- 
tional time scale, and it also showed radial dependence. 

An important dynamical feature of the time-dependent 
multi-arm spirals is that the spirals "co-rotate" with 
the galactic rotation. Figure |3] shows the angular pat- 
tern speed of a dominant stellar spiral (to = 4). The 
angular frequency is analyzed by the rate of change 
of phase 4> m =i{R) during two galactic rotation periods. 
The pattern speed is not constant. In fact, the pat- 
tern speed decreases with radius in a manner similar 
to galactic rotation for any radii. This result shows 
that there is no "single" co-rotating point in the disk, 
in contrast to stationary spiral waves (or bar) with 
a rigid rotation. Therefore, no single spiral can last 
for more than on e rot ational perio d (see Section [5]). 
IWada et all (|2011l ) and iGrand et all (|2012cQ have also 
reported that the spiral arms in their simulations are 
co-rotating, winding and short-lived. Observationally, 
the phenomenon of co-rotating spiral arms has been con- 
firmed in at least a few of the nearby spiral galaxies 
such as M51, NGC 1068, M101, IC 342, NGC 3938, and 
NG C 3344 using the Tremaine- Weinberg (TWR) analy- 
sis (jMerrifield et al.ll200fl iMeidt et al.ll200l l20M ). 

It is noteworthy that the pattern speed curve seems 
slightly flatter (particularly when R ~ 10 — 15 kpc) than 
the circular rotation curve. In order to investigate the 
cause of the flattening, we divided this period into five 
periods and performed the same analysis for each period 



(Figured]) ■ The pattern speeds in the all periods decrease 
with increasing radius in a manner similar to galactic 
rotation for any radii; however, they show slightly flatter 
distributions at certain instants (e.g., T rot = 12.4—12.8). 
This flatter distribution lasts for less than one galactic 
rotation period (AT rot < 1). Thus, this structure is not 
a long-lasting structure. 

Further, it is significant that the gravitational scatter- 
ing of stars by the spiral arms is not sufficiently large to 
erase all the spiral features. The curves in Figure [TJ; indi- 
cate that the value of Q increases from an initial value of 
1.2 to a final value of 1.4 around R « 2i? s d at T rot = 15. 
However, the spirals are not completel y erased due to dy- 
namical heating (see lFujii et al.l[2011| ). It is to be noted 
that heating effect along the vertical direction to the disk 
plane is negligibly small (bottom panels in Figure [2]), 
whic h is consistent with the results of previous stu dies 
(e.g.. lBinnev fc LacevlH988l : Uenkins fc Binnevlll990f ). 

4. GROWTH AND DAMPING PHASES OF STELLAR 
SPIRAL ARMS 

In this section, we investigate the dynamical evolu- 
tion of co-rotating spiral arms in detail. Figures [5] and 
[6] show the time evolution of a spiral arm in the co- 
rotating frame. A weak density enhancement can be ob- 
served around <j) ~ 240 - 300° and R ~ 7 - 10 kpc for 
Trot = 11-7. This density enhancement causes the growth 
of a prominent spiral arm until T rot — 12.0, and the arm 
has a maximum amplitude around T rot ~ 12.2 — 12.3. 
Beyond T rot ~ 12.3, this spiral arm rapidly fades out 
{T T ot = 12-3 — 12.5). The arm merges with with a neigh- 
boring weak spiral arm (T rot = 12.5 — 12.7), and the peak 
density contrast again becomes S ~ 1 when T rot = 12.8. 
In the following sections, we discuss the growth phase 
(Figure [5]) and the damping phase (Figure [5]) of the spi- 
ral separately. 

4.1. Growth Phase 

In the growth phase of the spiral arm (T rot ~ 12.0 — 
12.2), the pitch angle of the arm reduces due to the dif- 
ferential rotation of the disk. This is clearly observed in 
the right hand panels of Figures [5] and [6j the two peaks 
at 2i? s d + 1.5 kpc and 2i? s d — 1.5 kpc show increasing 
separation in terms of <j). More quantitatively, Figure [7] 
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Fig. 2. — (Top) Evolution of stellar disk. The disk's surface density is shown in the logarithmic scale. The rotation direction of the disk is 
counterclockwise. Here, a rotational period, T ro t, is measured at R = 2i? s( j, and corresponds to 325 Myr. (Middle) The radial distributions 
of the spiral modes analyzed by the Fourier decomposition (see text). (Bottom) The radial distributions of the vertical velocity dispersion 
(<t z ) of the stellar disk. 
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Fig. 3. — Pattern speeds of dominant mode (m = 4) for T ro t ~ 4.0 (left) and 12.2 (right). The pattern speed is analyzed by measuring 
the rate of change of phase (p m =4 at each annulus every 10 Myr, and these phases are combined over the periods T ro t = 3.0 — 4.0 and 
11.2 — 13.2. The contours show the amplitude, \A m —i{R)\, of the dominant mode. The superimposed curves show the radial variations in 

n,n± k/4. 
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Fig. 4. — Pattern speeds of dominant mode (m = 4) for T ro t ~ 12. The pattern speed is analyzed by measuring the rate of change of 
phase m= 4 at each annulus every 10 Myr during the periods T rot = 11.2 - 13.2, 11.2 - 11.6, 11.6 - 12.0, 12.0 - 12.4, 12.4 - 12.8, and 
12.8 — 13.2. The contours show the amplitude, \A m —4(R)\, of the dominant mode. The superimposed curves show the radial variations in 

n,n± k/4. 
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Fig. 5. — Evolution of spiral arm in growth phase. Left columns: face-on views of the disk. The surface density is shown in the 
logarithmic scale. The time interval for each phase of evolution is indicated in the upper corner of each panel to the left. The rotating 
frame at R = 2R sd = 8.6 kpc (indicated by a solid circle) is adopted as a reference. The dotted and dot-dashed circles indicate radii values 
of R = 2_R s( j — 1.5 kpc and 2i?. st j + 1.5 kpc, respectively. Middle columns: the density distribution corresponding to the spiral evolution 
along the (f> — R plane. Right columns: the corresponding azimuthal density contrast (S) profiles at R = 2_R s( j — 1.5 kpc (dot-dashed), 2iJ s( j 
(solid), and 2ii s( j + 1.5 kpc (dotted). 
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Fig. 6. — Evolution showing thinning out of spiral (T ro t = 12.2 — 12.5), with subsequent rcconncction with neighboring spiral arm 
(T rot = 12.5 - 12.7). 



8 



Baba, Saitoh & Wada 




Arm Pitch Angle i [deg,] 



Fig. 7. — Evolution of spiral arm on i — <5 plane for T ro t = 12.0 — 
12.5. The hatched region corresponds to the predicted maximum 
pitch angle around the analyzed region (Q ~ 1.4 and T ~ 0.75 — 
0.85) due to swing amplification (refer to equation (7)). 



shows the evolution of the spiral arm along the pitch an- 
gle - density contrast (i — S) plane. Here, we calculated 
the pitch angle (i) of the spiral arms using the relation 



1 AR 

tarn = — . 

R A0' 



(6) 



where R = 2i? sc i (= 8.6 kpc), AR = 3 kpc, and A(f> is the 
azimuthal angle difference between the contrast peaks at 
R = 2R s d — 1.5 kpc (dot-dashed lines) and R = 2i? s( j + 
1.5 kpc (dotted lines) shown in Figures [5] and [6] Here, 
we evaluated the positions of the contrast peaks at each 
radius by visual inspection. The arm density contrast, 
5, is calculated by averaging the corresponding density 
contrast over the radial range in Figures [5] and [51 As 
the pitch angle of the spiral arm decreases from i ~ 40° 
(Trot = 12.0) to i ~ 32° (T rot = 12.20), the density 
contrast increases to a maximum, and subsequently, it 
decreases with increase in the pitch angle. Thus, the 
spiral arm has a maximum amplitude when i ~ 32°. 

This amplification process associated with the galac- 
tic she ar motion is known as the swing amplifi- 
cation dToomTi 119811; jGoldreich fc Lvnden-Belll 119651: 



Julian fcTbomrd 119661: iToomre fc Kalnaisl 11991 " 
Fuchsl (|2001f ) calculated grid models of swing amplifica- 
tion process by varying the shear rate T for Q = 1.4. In 
his work, he derived a fitting empirical equation for swing 
amplification (eq.98 in his paper) as the following: 



tanith = 1-932 - 5.186 



4.704 



(7) 



where ith is the pitch angle at which the spiral arm 
reaches the maximum amplitude. We adopt this equa- 
tion to evaluate the predicted pitch angle i t h- In our 
model, T ~ 0.8 around R = 2R sd (Figure Thus, 
substituting T ~ 0.8 in eq. (7), we obtain i t h — 32°. 



5 D'Onghia ct al. (2012) have investigated the growth of spiral 
arms via swing-amplification, and their nonlinear evolution is not 
ful ly consistent with the cla ssic swing-amplification picture proved 
by I Julian fc Toomrel [Tl966h . 



This value is consistent with the evolution of the spi- 
ral arm with the maximum contrast (hatched region in 
Figure [7]). 

4.2. Damping Phase 

The swing amplification mechanism can explain cer- 
tain aspects of the evolution of spiral arms, i.e., the am- 
plification (excitation) of density enhancement. How- 
ever, the destruction proc ess of the non-steady spirals, 
as seen in Figure IS (see also lFujii et aD20ntlWada et alJ 
1201 lh . cannot be understood only by the swing amplifica- 
tion mechanism. The top right panels in Figure [5] show 
the time evolution of the relative velocities (Av^ and 
Av^) in the co-rotating frame of the spiral arm. Upon 
considering the coordinates (?7,£), i.e., the r\- and £-axes 
are perpendicular and parallel to the spiral arm, respec- 
tively, the velocity component along each is given by 

Aw,, = vr cosi + (vf - V c i r ) sim, (8) 

Av£ = —v R sin % +(v^ — v C ir) cos i, (9) 

where i is the pitch angle of the spiral arm, and w c i r is 
the circular velocity determined by the azimuthally aver- 
aged gravitational field. For this definition, inflow to the 
arm corresponds to Av v > for rj < 0, and Av v < for 
T) > 0. In the amplification phase (T rot < 12.2), we see 
a clear inflow motion to the spiral arm along both sides 
of the chosen strip (solid lines). In the initial stages of 
the damping phase (T rot ~ 12.3), the streaming velocity 
around the spiral arm gradually transits from inflow to 
outflow. At the end of the damping phase (T rot ~ 12.5), 
it is clear that the stars in the spiral arm move away from 
the arm along both its sides (dot-dashed lines). During 
the dynamical evolution, the parallel component of the 
velocity, Av^, does not change its sign or show a de- 
crease its magnitude. The left panels of Figure [5] show a 
schematic view of the time evolution of the non-circular 
velocity field associated with the non-steady spiral arm. 

The middle right panels in Figure [5] show the time 
evolution of the non-axisymmctric gravitational forces 
(i.e., spiral perturbation) involved in the damping phase. 
The component of the non-axisymmetric gravitational 
force perpendicular to the spiral arm (A-Fg raVj ,,), which 
is stronger than the parallel component (A.Fg raVi £), is al- 
ways directed towards the arm (A_Fg rav ^ > for r\ < 0, 
and AFgrav.jj < for i] > 0) during the dynamical evolu- 
tion. The bottom right panels in Figure [H] show the time 
evolution of the net force (i.e., non-axisymmctric gravita- 
tional force plus Coriolis force perturbation). The com- 
ponent of the net force perpendicular to the spiral arm 
(AFg rav _|. C ori ,-q) evolves in the same manner as AF graVit7 
because there is a strong density gradient along the rj- 
direction. The parallel component A_Fg rav ^ is almost 
zero during the growth and damping phases and does not 
change. This is because the density gradient along the 
^-direction is small and almost unchanged. However, the 
parallel component of the net force Ai^grav+cori^ changes 
its 77-dependence from the growing phase (T rot = 12.1) to 
the damping phase (T ro t > 12.3) according to the change 
of the sign of the perpendicular velocity (Ai^), since the 
Coriolis force works perpendicular to the direction of the 
velocity. 

This indicates that the Coriolis force exerted on the 
stars in the damping phase exceeds the non-axisymmetric 
gravitational force due to the spiral perturbation. This 
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Fig. 8. — (Left) Schematic view of velocity field with respect to non-steady spiral arm. The spiral arm is represented by the gray region. 
(Right) Top row: time evolution of relative velocities with respect to the stellar spiral arm (~ 270°), At),, (left panel), Avg (right panel), at 
T ro t = 12.1 (solid line), 12.3 (dashed line), and 12.5 (dot-dashed line). The analysis is performed for the strip along the »7-axis with a width 
of |£| = 500 pc. Middle row: time evolution of non-axisymmetric gravitational forces (i.e., spiral perturbation) AF gTlLV = F gIav — -F c i r . 
Here, F clI denotes the axisymmetric gravitational force calculated by the azimuthally averaged gravitational forces. Bottom row: time 
evolution of "net" non-axisymmetric forces AF gmv+cor i(= AF gT!iv + Ai^o,-;). Here, the Coriolis force perturbation AF COI i is calculated as 
— 2S7 X (v — 1> c ir). 



causes stars to "escape" from the spiral perturbation, 
and eventually the spiral arm itself begins to thin out 
and fade. 

The above argument suggests that the non-steady na- 
ture of stellar spirals is originated in the evolution of the 
orbits of the stars i n the spirals. The phenomenon of 
swing amplification (|Toomrelll981[) is a part of this non- 
linear coupling between particles and waves, but it does 
not describe all observed phenomena. In the next sec- 
tion, we explore the orbital evolution of stars associated 
with the spiral arm in detail. 

5. ORBITAL EVOLUTION OF STARS AROUND SPIRAL 
ARMS 

Figure [5] shows the evolution of stars along the <f> — R 
plane in the early phase of spiral development. At this 
stage of the simulation, we selected 15 particles associ- 
ated with one of the three weak spiral arms that had 
evolved at T TOt = 4.0. The figure shows that stars with 
epicycle motion arc captured by the density enhancement 
(Trot = 3.6 — 4.0), and further, these stars are dissociated 
from the original arm (T rot = 4.2 — 4.6). This behavior is 
similar to that of the "density wave" from a certain view- 
point, but the spiral arms are short-lived and co-rotating 
(left panel of Figure [3]). Thus, even if the arms appeared 
in the early linear phase, they are not completely ex- 



plained by the picture of stationary density waves. The 
middle column of images in Figure [S] shows the plots of 
the azimuth angle {<f>) versus the angular momentum L z 
curve instead of the <f> — R plot. It is clear that all the 
stars along the <j> — L z plane oscillate horizontally when 
the angular momentum of each star is conserved. The 
images in the right column of Figure [9] show the so-called 
Lindblad diagram, where the angular momentum L z of 
each star is plotted against its total energy E. Most of 
the stars show no significant movement from their origi- 
nal position, thereby suggesting that their energy is also 
conserved. 

On the other hand, the behavior of the stars is very 
different in the non-linear phase, where the spiral arms 
arc well-developed and non-steady (refer to previous sec- 
tions). Figure [TU] shows images identical to those in Fig- 
ure [9j however, in this set of images, the stars are in 
motion due to change in their angular momenta in the 
non-linear phase of orbital evolution (see also the sup- 
plementary video). When the stars are captured by the 
density enhancement (T ro t — H-8 — 12.0), they radially 
migrate along the spiral arms. The stars approaching 
from behind the spiral arm (i.e., inner radius) tend to 
attain increased angular momenta via acceleration along 
the spiral arm, whereby they move to the disk's outer 
radius. In contrast, the stars approaching ahead of the 
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spiral arm (i.e., outer radius) tend to lose their angular 
momenta via deceleration along the spiral arm, and they 
move to the disk's inner radius. Along the <$> — L z plane, 
the stars oscillate both horizontally as well as vertically. 
Moreover, the guiding centers of the oscillations do not 
remain constant at the same value of L z . This is essen- 
tially different from the epicycle motion in which L z is 
conserved. 

The Lindblad diagram in Figure QJJ shows that the 
stars oscillate along the curve of circular motion by un- 
dergoing change in terms of both angular momentum and 
energy. The oscillating stars successively undergo aggre- 
gation and disaggregation along the curve, thereby lead- 
ing to the formation of structures referred to as "swarms 
of stars" along the <f> — L z and R — cf> planes. The non- 
steady nature of the spiral arms originates in the dynam- 
ical interaction between these "swarming" stars with a 
non-linear epicycle motion and the high-density regions, 
i.e., the spiral arms moving with the galactic rotation. 
This is similar to the wave-particle interaction described 
previously; however, in this case, the high-density regions 
are not "waves". Since the non-steady spirals move at 
the rate of the local galactic rotational speed (Section 
3), in contrast to a spiral perturbation with a single pat- 
tern speed, co-rotating points are found everywhere on 
the spiral arms. Therefore, the motion of stars along 
the E — L z plane can be naturally understood as due to 
their scattering around co-rotating spirals 6 . This be- 
havior is similar to reported in previous studies (e.g., 
Sellwood fe Binneyll20"0llGrand et aTll2012d : lBi"rd et al.1 
2012tlGrand et al.ll2012al) . Furthermore, we observe that 
the E — L z curve of each star changes over a large range 
of radii. This is entirely different from what is expected 
in stationary density waves, whe re these changes are lim- 
ited t o the Lindblad resonances (|Lvnden-Bell fc Kalnaisl 
11972ft . Further, it is noteworthy that the structures 
formed self-induced in the angular momentum space; this 
is a property similar to that of the "groove" mode hy- 
pothesized by ISellwood fc Linl (|1989D . This point is be- 
yond the scope of this paper. We will investigate this 
elsewhere. 

The time evolutions of the angular momentum of stars 
associated with a spiral arm at T rot ~ 4.0 and T rot ~ 12.2 
are plotted in the left top panels of Figures [TTJ and fT2l 
respectively. In the early phase (T rot ~ 4.0), both the 
angular momentum and random energy do not change by 
more than 10% during one rotational period (top right 
panel of Figure [TT]) . On the other hand, in the non- 
linear phase (T rot ~ 12.2), the fraction of the angular 
momentum changes by ~ 50% (top right panel of Figure 
[T2|) . It is clear that the angular momentum of the stars 
changes significantly due to their scattering by a well- 
developed spiral at T rot = 12.2. This corresponds to the 
hypothesis that the guiding center of the epicycle motion 
of each star undergoes radial motion. 

The change in the normalized random energy shown 
in the left bottom panels of Figures [TTJ and [12] indicates 
that certain stars with relatively small initial random 
energies experience a large energy change Ai? ran( j after 
interaction with a spiral arm. Here, the random energy 

6 Stars around the co-rotation point change their an- 
gular momenta without increasing their random energy 
ifLvnden-Bell fe KalnalillT972T ) . 



is calculated as the difference between the total energy 
and the circular energy, i.e., E vaxi & = E — E cu (R gc ). It is 
noteworthy that the random energy change is not always 
positive; a significant fraction of stars lose their random 
energy. This is because the perturbation from the spiral 
arm shifts the guiding center of the stars' epicycle motion 
without increase in orbital eccentricity. It is to be noted 
that this argument is rigorously correct; in fact, Figure 
[T"2l (bottom left panel) shows the changes in the random 
energy of the stars as a function of change in the angular 
momentum. We can see a weak trend: the outer migra- 
tors (i.e., ifin — L- m i > 0) undergo a decrease in their 
random energy (.E ra nd,fin - £rand,ini < 0), and vice versa. 
A similar effect of the radial migration of stars around 
spiral arms upon disk heating has been noticed in recent 
numerical simulations ([Grand et al.ll2012d | Rqskar et al.l 
[20TllGrand et al.ll27jl~2a! iMinchev et al. f l2012fl . 

In summary, the gravitational interaction between the 
stars in the spiral arm and the spiral density enhance- 
ment changes the angular momentum and random en- 
ergy of the stars, and this process in turn changes the 
structure of the spirals. During this process, the ran- 
dom energy of individual stars in the system does not 
increase monotonically. In other words, local interac- 
tions between the non-steady arms and stars increase 
or decrease the total energy of individual stars locally; 
however, the energy remains around its value for circular 
motion with the occurrence of a small dispersion. This 
is because the interaction causes the migration of the 
guiding centers of the stars without increasing their ec- 
centricity or random energy. This "dynamical cooling" 
mechanism is essential to prevent heating of the stellar 
disk and erasure of the spiral arms, and the mechanism 
produces "swarms" of stars moving between non-steady 
spirals. The non-linear epicycle motion of the stars and 
their non-linear coupling with the density perturbation is 
the fundamental physics of the recurrently formed, non- 
steady spiral arms in a stellar disk. 

6. DISCUSSION 

Figure [Jj2 compares the morphological properties (the 
pitch angle i and the amplitude |^4m|) of the stellar spi- 
ral arms between simulations and observations. For the 
comparison, we used a 2D fast Fourier transform (2D 
FFT) a nalysis with logarithmic spirals (for details, see 
for e.g.. iPuerari fc Dottoril[l99 2[). The plotted points in 
green and red symbols represent the time evolution of the 
spiral arms for the two models that are essentially iden- 
tical (see footnote). We found that the simulated spiral 
arms exhibit a maximum amplitude around i ~ 30° . Fur- 
ther, we found that the spiral arms tend to be weaker in 
the pure stellar disk model (model MS) than the model 
with the ISM (model MSG) 7 . This distribution which 
has a peak at around i ~ 20 — 30° is qualitatively con- 
sistent with the the prediction of the swing amplification 
mechanism (see section 4.1). We compared the simulated 

7 Mode l MSG is the same model as that presented by 
IWada c t al. (201]]). This model is mostly identical to model MS; 
however, in model MSG, the ISM, star formation, and supernova 
feedback are also taken into account. The hydrodynamics of the 
model is solved by the smoothed particle hydrodynamics (SPH) 
method. The initial gas mass fraction is 10 % of the stellar disk 
mass, and the initial density profile follows an exponential pro- 
file with a s cale-l ength twice that of the stellar disk. Refer to 
IWada et all l20lTI ) for details. 
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Fig. 9. — Orbital evolution of stars in spiral arm. The stars associate around the spiral arm within a distance of ±0.5 kpc at T ro t = 4.0. 
Left columns: orbits on <f> — R plane. Middle columns: orbits on <p — L z plane. Right columns: orbits on E — L z plane. The colors denote 
the angular momentum at the time instants when the stars associated with the spiral arm. See the attached video. 
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Fig. 10. — Orbital evolution of stars in spiral arm during non-linear phase. The stars associate with the spiral arm at T ro t = 12.2. See 
the attached video. 
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Fig. 11. — Time evolution of orbital characters of stars around spiral arm during early evolutional stage shown in Figurc[9] Left: Time 
evolution of the orbital angular momentum (top) and the random energy (bottom) of the stars. The plotted stars are the same as those in 
Figure [9] The random energy of the stars is non-dimensional using the radius of the guiding center R gc and epicyclic frequency K gc , which 
is proportional to the square of the orbital ellipticity e. Each star is indicated by a different line thickness. The colors correspond to those of 
Figure |91 Right: Angular momentum change, Lfl n — L- ln i, distributions of stars (top). Here, AT ro t = 0.8 is adopted to measure the angular 
momentum change and the vertical axis indicates the normalized number of particles. The random energy changes £? r and,fin — ^"rand inii m 
the stars over AT ro t = 0.8 as a function of their angular momentum changes (bottom). The colors correspond to those in the left panels. 



distribution with the observed one. The dependence of 
\A m \ on the pitch angle in the models with/ without the 
ISM is consistent with observations at least for the range 
of values corresponding to i < 30°. This comparison sug- 
gests that the self-induced spiral arms in differentially 
rotating disks and their time evolution can fairly con- 
sistently account for the morphological diversity seen in 
nearby spiral galaxies. 

However, this comparison between the simulated 
and observed distributions requires further examination. 
Firstly, the determination of the pitch angle is some- 
what uncertain. The analyzed pitch angles differ be- 
tween studies for the same galaxy. For example, the 
pitch a ngles of NGC 3054 ar e given by 3 3 °, 43° , and 
12° in lGrosb0l et all (l200l . iSeigar et al.l (l2006lh and 
iDavis et al.l ( 2012|) . respectively. Secondly, there is a lack 
of statistical studies on the relation between the pitch 
angle and the amplitude of the stellar spiral arms. Al- 



though some statistical studies have used 2D FFT to 
analyze the pitch angle of the spiral arms, they do not 
include examinati on of the amplitud e or contrast of the 
spira l arms (e.g., ISeigar et al.fl2005l 120061; IDavis et al.l 
120121 ) . Thus, further progress in understanding the spi- 
ral dynamics requires additional statistical and robust 
observation data of the morphological properties (both 
the pitch angle and amplitude) of the stellar spiral arms. 
This relation between the pitch angle and amplitude will 
form one of the tests for spiral genesis theories along 
with the pitch angle-shear rate (or Hubble-type) relation 
dRoberts et al| l975: Sei gar fc Jamesll998HHozumill2003l : 



iGrand et al.ll2012bl ). existence of systematic angular off- 
set between the young stellar component and the stellar 
spira l arm (iFuiimotol 1968t iRoberti |1969t lEgusa et all 
l2009HFoyle et al.ll201lt iFerreras et al.ll2012D. and the ra- 
dial d ependence of the pattern speed (jMeidt et al.ll2008l 
2009). 
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Fig. 12 . — Time evolution of orbital characteristics of st ars f or phase in which the spiral arms are non-linearly developed (shown in 
Figure [TOt . The plotted stars are the same as those in Figure ITDl 



Finally, we co mment on grand - design spirals. In a se- 
ries of studies (|Fuiii et all 1201 H IWada et all 120111 ) in- 
cluding this one, we have examined the stellar dynam- 
ics of non-steady stellar spirals, as well as the interac- 
tions between them and the ISM. However, these studies 
have focused more on multi-armed spirals than grand- 
design spirals (i.e., m = 2 spirals) whose fr action is more 
than ~ 50% in nearby spiral galaxies (|Grosb0l et all 
120041 Kendall et all 120111 ). It has been observed that 
gran d-design spirals are associated with bars or compan- 
ions (iKormendy fc Normanl 119791: ISeigar k Jameslll998l 
iSalo et al.ll2010t Kendall et all 120111 ), and this observa- 
tion is consistent with the results of many numerical sim- 
ulations of bar - driven spirals (e.g., ISellwood k Sparki 



19881: iBotlemal 
Oh et all l2008t 



20031) and t idally - driven spiral s (e.g 



Dobbs et all Eoift LStruck et all Eoill ). 



We will focus on bar-driven spirals and tidally-induced 
spirals in forthcoming studies. 

7. SUMMARY 

The TV-body simulations of an isolated disk galaxy 
show the formation of self-induced, non-steady multi-arm 



spirals that follow the differential galactic rotation. Wc 
found that the swing amplification mechanism causes the 
development of the spirals. When a spiral undergoes the 
damping phase, the Coriolis force dominates the grav- 
itational perturbation exerted by the spiral, and as a 
result, stars escape from the spirals, and join a new spi- 
ral at a different position. This process uniformly for 
a given spiral, thereby resulting in the formation of bi- 
furcating and merging spiral arms; therefore, the domi- 
nant spiral modes always show change in their radii over 
time. We confirmed that this phenomenon originates due 
to the changing orbital properties of stars. The angu- 
lar momentum and energy of each star undergo changes 
due to the star's interaction with the spiral arms. As 
a result, the epicycling stars radially migrate; in other 
words, their guiding centers also undergo motion. Inter- 
estingly, the movement of groups of stars with similar 
orbital properties causes the appearance of "swarming" . 
In the non-linear phase of development of spiral instabil- 
ity, the swarming stars cause complicated morphological 
changes in the spiral arms. 
During this process, the random energy of individual 
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Fig. 13. — Pitch angle-amplitude correlation. Morphological 
properties (amplitudes A m and pitch angles i) of stellar spiral arms 
are derived using a 2D Fourier fitting to the range of radii given 
by R a d < R < 3R s d. Here, \A m \ values for the m = 4 mode in the 
simulations are shown. Model MSG is nearly identical to Model 
MS; however, model MSG also includes the gaseous component 
with a mass fract ion that is 10% of the stellar disk mass (refer to 
IWada e t al. 2011!). The filled tri angles represen t the o bservational 
data obtained from the study by Grosb0l et al. (2001). 



stars (or orbital eccentricity) does not increase monotoni- 
cally. In fact, a significant fraction of stars even lose their 
random energy. This "dynamical cooling" due to the 



mechanism like the wave-particle interaction can explain 
why the short-lived spiral arms are self-induced over sev- 
eral rotational periods despite the absen ce of dissipative 
component in the disk ([Fuiii et al.ll2011lh 

In the above process, it is essential that spiral arms 
mostly follow the galactic rotation at any radius; in 
other words, the "co-rotating points" are required to 
be ubiquitous in the differentially rotating galactic disk. 
The conclusions of other previous studies also indi- 
cate the possibility that the co-rotation resonance af- 
fects stellar motions more in terms of radial migra- 
tion than heating up of the disk (Lvnden-Bell fc Kalnai 
1971 ISeilwqod fc ; Binnevl l200l HoTand et al.l l2oT2cF 
Roskar et alj|2012t iMinchev et al.ll2012l ). 

We conclude that the non-linear epicycle motions and 
self-gravity in the differential rotation of stellar disks are 
essential for the recurrent amplification and destruction 
processes of the spiral arms. In other words, the issue of 
the so-called "winding dilemma" is no longer a problem 
at least in multi-armed spiral galaxies. 
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